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ABSTRACT 


Narrow-waveband (lOOA) photoelectric slit scan photometry of the 
Neptune disk is reported^ Observations were concentrated within the 
strong CH/, band at X7300A. For comparison, measurements were also made 
within a continuum waveband at A6800A. Point spread function data were 
obtained in both colors. Qualitative estimates of the true intensity 
distribution over the Neptune disk were made. Within the X6800A continuum 
band, Neptune appears as an essentially uniform disk. Within the X7300A 
CH^ band, the planet exhibits strong limb-brightening. Our results appear to 
require the presence of an optically thin layer of brightly scattering 
aerosol particles high in the Neptune atmosphere. 
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1 . INTRODUCTION 

Trafton (1974) reported the first measurements of the H 2 quadrupole 
lines in the spectrum of Neptume. Strengths of the S(0) and SQ) lines of 
the (4-0) band were interpreted in terms of both an inhomogeneous atmosphere 
overlying a reflecting layer, and a homogeneous, semi-infinite scattering 
atmosphere. Only the scattering model proved to be consistent with Neptune's 
spectrum in this wavelength region. The H 2 abundance along a scattering 
mean free path was found to be 'u 320 km. amagat; the maximum permitted 
abundance was 450 km. amagat. By comparison, the scattering mean free 
path for Rayleigh .-scattering in a pure H 2 atmosphere is equivalent to 715 
km. amagat. Although aerosol particle scattering may be responsible for 
the observed reduction in the mean free path, Trafton noted that an al- 
ternative explanation might be molecular scattering from gases additional 
to- hydrogen. From an analysis of photoelectric photometry of Neptune, 
Lockwood C1978) concluded that the scattering properties of the atmosphere 
are variable, and that aerosol particles are present. 

Infrared (1 -4 pm) observations of Neptune have been interpreted 
by Joyce et al . (1977) and by Pilcher (1977) to suggest that, over a one 
year period, an extensive high altitude cloud formed on the planet and 
then partially dissipated. About 10-months after its initial appearance, 
the optical thickness of the cloud averaged over the planet was on the 
order of unity. An H 2 column-abundance of 5+2 km. amagat was derived; 
it appeared to refer to an effective reflecting layer high in the atmosphere 
of Neptune. 

Recent infrared (0,8 - 2,5 pm) studies of Neptune by Fink and 
Larson (1979) have indicated that its atmosphere has a variable height 
cloud layer overlain by 50 km. amagat of H 2 . Fink and Larson point out 
that the high cloud layer cannot have a large optical thickness since both 
the H 2 quadrupole lines and the visible CH^ bands give quite large abundances. 
Specifically, the CH^ column abundance was found to be 0.7 km. amagat. 

For Uranus, the corresponding value was 1.6 km. amagat. Fink and Larson 
also found the Uranus atmosphere to be clear to great depths with its upper 
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regions depleted in CH^. It appears likely that visible radiation can 
probe to great depths in the Neptune atmosphere. 

Useful information concerning the opacity of the Neptune 
atmosphere, in particular the presence of an optically thin upper cloud 
layer, can be derived from measurements of the distribution of intensity 
over the disk of the planet both within and between the strong CH^ bands. 

In this paper, we report the first attempt to infer the wavelength-dependent 
optical appearance of the Neptune disk. 


i 


2. OBSERVATIONS 

0 

During the 1979 Neptune apparition, narrow-band (lOOA) slit scan 
photometry of the planet was carried out within two regions of the spectrum 
. by the technique described by Price and Franz (1979). Scans were obtained 
primarily within the 730o8 CH^ band. But, for comparison, scans were also 
made within the continuum region at 680o8. Waveband selection was achieved 
by means of two filters used in our early investigation of Uranus (Price 
and Franz: 1976). On 1979 July 23, reliable observational data were obtained 

with the area-scanner mounted on the 1.8 meter Perkins reflector at Lowell 
Observatory. 

Neptune was difficult to observe because of its extreme southerly 

declination C-21 degrees). Even at meridian transit, the planet never 

reached an elevation angle greater than 'v 33 degrees. In order to minimize 

’the path length through the atmosphere, observations were restricted to 
* * 

within ± 2 hours of transit. Effects of differential refraction, both 
within and between each waveband, were avoided by orienting the slit length 
in a, north-south direction. Both the stellar and Neptune images would then 
suffer dispersion along the slit when observations were made near the 
meridian. Special care was taken in the determination of the point spread 
function. A star, of apparent visual magnitude 'V' 7.0, located within 
45 arc. min. of Neptune, at an elevation angle essentially identical with 
that of the planet, was selected to provide PSF information. Both Neptune 
and the star were scanned in an east-west direction only. Both forward 
and reverse scans were used. Samples were taken at spatial intervals of 
0.032 arc. sec with a slit of 100 -pm (0.645 arc. sec) width. Both slit and 
scan length were 2mm (12.9 arc. sec). 

Each recorded scan of Neptune and of the star was obtained by 
integrating 100 one-second scans. Shortening the integration time, although 
desirable to minimize image wander produced by telescope tracking errors 
and by variations in atmospheric refraction, would have adversely affected 
the signal-to-noise ratio. If image wander did occur, we hoped to match 
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its effects for both Neptune and the star by equalizing their integration 
times. To increase the effective signal-to-noise ratio, composites of the 
integrated scans were subsequently computed by the method described by 
Price and Franz (1979). Scans made in opposite directions were mirrored 
on their centroids prior to computer summation. Both the Neptune and 
stellar composite scans were symmetrical within the level set by the 
residual photometric noise. Total scan numbers for both the Neptune and 
stellar composites are shown in Table I. Final Neptune and stellar slit 
scans were formed by averaging the respective composites about their centroids. 

Atmospheric turbulence, together with diffuse instrumental scattering, 
•is responsible for the observed point spread function. Our earlier investi- 
.gations of Uranus (Price and Franz: 1978, 1979) showed that the point 

spread function is well represented by the summation of two colocated 
but distinct Gaussian curves. Normalizing the point spread function, f(r), 
to unit total energy* we can write its radial distribution in the form 

f(r) = 2 ^ 2-, + B exp (1) 

ir [Aa^ -f Bog ] ( ' * ) 

where A, B, and are constants. Optimum sets of PSF parameters for 

both wavebands were obtained from the composite stellar slit^^cans by the 

method described by Price and Franz (1979). Within the 6800A v.'aveband, 

we found A = 1, B = 0.1336, a, = 0.9829 arc-r sec, and = 2.0730 arc, sec. 

Within the 7300A waveband, we found A = 1, B = 0.1378, = 0.9829 arc. sec, 

and - 2.0413 arc. sec. As one would expect, the point spread function 

became narrower with increasing wavelength. Our best-fitting analytical 

scans are compared with the observational data in Figures 1 and 2 for the 
0 o 

6800A and 7300A wavebands, respectively. Scaled with respect to each 
central intensity, the respective root mean square deviations between 
calculation and observation are 0.62 percent and 0.68 percent. Evidently, 
the goodness-of-fit in each case is limited by the residual noise in the 
photometric data. Stability of the point spread function during the 
course of the observations was high. Consider both sets of PSF data 
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Table I 


FORMATION OF COMPOSITE STELLAR AND NEPTUNE SLIT SCANS 


Waveband 

(Angstroms) 


Integrated 

Scans 

One-Second 

Sweeps 

6800 

star 

2 

200 

Neptune 

4 

400 

. 7300 

Star 

2 

200 

Neptune 

8 

800 


8 


















together. The Gaussian 1/e-widths of the individual integrated stellar 
slit scans exhibited an rms deviation equivalent to 2.0 percent of the 
mean. Residual uncertainty in the mean point spread function for each 
waveband will have a negligible effect on the interpretation of the 
planetary data. 
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3.. ANALYSIS 


Careful examination of the Neptune slit scan data showed the level 
of the residual photometric noise to be too high to permit reliable restora- 
tion of the actual radial intensity distribution over the disk. Instead, 
we chose to qualitatively Investigate the appearance of the planet by 
comparing the observed slit scans directly with predictions made for a 
variety of model intensity distributions. Our approach was first to model 
both the size and shape of the planet and to adopt a candidate radial dis- 
tribution of intensity over the disk, next to employ the known point spread 
function in a two-dimensional convolution to derive the planetary image 
smeared by atmospheric seeing, then to compute the profile which would 
result from slit scanning the image in one-dimension (E-W), and finally 
to normalize the slit scan prediction to permit comparison with the observa- 
tional data. Mathematical procedures developed by Price and. Franz (1978) 
were used. 

From stellar occultation data. Freeman and Lynga (1970) found 

Neptune to be oblate, with an equatorial radius of 25265 +36 kms and an 

ellipticity of 0.026 +_ 0.005. Their data refer to a level in the atmosphere 

corresponding to a molecular number density in the range 10 to 10 cm . 

Assuming the Neptune atmosphere to be isothermal, Freeman and Lynga derived 

an equatorial radius of 24753 + 59 kms at a molecular number density of 
20 -3 

10 cm where they expected a dense cloud layer to form. But, whether 
or not the Neptune atmosphere is indeed clear to such a depth is uncertain. 
For our present analysis, we adopted an equatorial radius of 25,000 kms. 
Effects of oblateness remain to be evaluated. Detailed consideration 
of both the observational geometry and technique showed that accurate 
theoretical slit scans of Neptune could be derived if the planetary 
disk were assumed to be perfectly circular with a radius equivalent to 
0.994 times the equatorial value. Since the Neptune radius, relevant 
to our analysis, is not known within an accuracy of 0.5 percent, we 
chose to ignore the correction for oblateness. Neptune was taken to 
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be perfectly spherical with a radius of ?,5,000 kms. The American Ephemeris 
and Nautical Almanac gives the geocentric distance of Neptune on 1979 duly 23 
as 29.515 AU. The corresponding angular radius of the planet was 1.168 
arc. sec. 

Six simple models were adopted for the actual distribution of 
intensity over the Neptune disk. Circular symmetry was adopted throughout. 

In five models, the intensity was assumed to be a linear function of 
radial distance from the center of the disk. Each model was defined by the 
ratio of the intensities at the center and limb of the planet (CTL). A 
complete range of intensity distributions, from strong limb darkening 
through a flat disk to strong limb brightening, was considered. A sixth 
model represented the most extreme sicuaticn nf limb brightening with the 
intensity assumed to be zero er-jept within a ring of infinitesimal width 
located at the limb. Characteristics of the six models are summarized in 
Table II. 

Theoretical Neptune slit scans are compared 1n Figure 3. All six 
models have been normalized to the same integrated energy. In general, the 
computations were based on mathematical procedures developed by Price and 
Franz (1978). But, for the ring of infinitesimal width, the radial intensity 
distribution of the smeared image was calculated by the method described in 
an Appendix. For both wavebands observed, the point spread functions are 
very similar. We therefore adopted a mean PSF for the comparison of the 
individual models. Specifically^ its parameters were A = 1, B = 0.1357, 
aij = 0.9829 arc. sec and 02 = 2.0572 arc. sec. In spite of the small 
angular diameter of the Neptune disk, strong limb-darkening or limb- 
brightening would be detectable if present. Note that the predicted scan 
intensity is essentially independent of the disk model at Ml. 8 arc. sec 
from the profile center. 

In Figures 4 and 5, observed^Neptune slit scan data are compared with 
theoretical predictions for the 6800A and 7300A wavebands respectively. 

To avoid cluttering the diagrams, theoretical predictions are plotted 
only for the flat disk and for the ring of infinitesimal width. Predic-: 
tions were based on the individual point spread functions relevant to each 
waveband. Both the observational data and theoretical predictions have 
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Table II 


MODEL INTENSITY DISTRIBUTIONS FOR THE NEPTUNE DISK 


Model 

CTL 

1 

1;0 

2 

1 r 0 . ■.> 

3 

1;1 

4 

0.5:1 

5 

0:1 

6 

Ring 
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been normalized to the same integrated energy. By inspection, we conclude 
that Neptune is well represented by flat disk at 6800A, while strong 
limb-brightening is clearly present v/ithin the 7300A CH^ band. The degree 
of limb-brightening, however, cannot be determined accurately because of 
residual photometric noise. 

Where the profiles are insensitive to the disk model, the observed 
and theoretical scans are very well matched. This consistency is most 
significant. It implies that tracking errors were negligible during 
acquisition of the Neptune data, that the radius of the planet was properly 
chosen, and that the observed point spread function provided an accurate 
description, of atmospheric smearing of the disk. 
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4. DISCUSSION 


Limb-brightening on Neptune must be the result of a vertical 
inhomogeneity in its atmosphere. More specifically, the single scattering 
albedo must decrease inversely with optical depth. Within the 730QA 
CH^ band, two explanations of the limb-brightening phenomenon are feasible. 

In the absence of aerosol particles, limb-brightening may be the result of 
methane depletion in the upper regions of the atmosphere. But, if the 
CH^/H 2 mixing ratio is constant with optical depth, limb-brightening may 
be caused by an optically thin layer of aerosol particles high in the Neptune 

atmosphere. HybHd explanations are also possible. 

0 * 

Within the 6800A continuum waveband, limb-darkening would be ex- 
pected to occur if the atmosphere were homogeneous. In a clear Hg - CH^ 
atmosphere, vibrational Raman scattering from molecular hydrogen would 

reduce the monochromatic single scattering albedo significantly below 

0 

unity. ’Absorption from a weak CH^ band at 6825A would also contribute 
to ii reduction in the single scattering albedo. Within this methane band, 
Giver (1978) found pressure effects to be insignificant. So no vertical 
inhomogeneity would result. For Neptune to appear as an essentially 
uniform disk, a limb-brightening mechanism must be acting to counter the 
natural tendency to limb-darkening. Such a mechanism could well be the 
presence of an optically thin layer of brightly scattering aerosol particles 
at high altitude. In fact, to reconcile our observational results with 
those obtained by other techniques, such a scattering layer appears to be 
essential. But whether or not the aerosol haze is laterally homogeneous 
cannot yet be determined. 
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APPENDI.X 

ATMOSPHERIC SMEARING OF A RING OF INFINITESIMAL WIDTH. 


Consider the situation of extreme limb brightening on a perfectly 
circular planetary disk. Suppose that all the radiation reflected By 
the planet is concentrated within a ring of infinitesimal width and with 
a radius equal to that of the disk itself. Normalizing the total radiant 
energy to unity, we can define the radial intensity distribution, g(r), 
in terms of the delta function. Specifically, we can write 

• ^ 

where a is the radius of the’ planetary disk. Suppose further, that atmospheric 
turbulence produces a smear described By a circularly symmetric point 
spread function, f(.r)» of the form 

f(r) — J .{a exp [-r^/o,^J + B exp [-r^/op^j ! (A-2) 

ir [Aa^ + Bap j ( • ' J 

where A, 6, and are constants which may be determined from observational 
data. Convolution of the point spread function, f(r), with the degenerate 
annulus, g(r), will produce a circularly symmetric image whose radial 
intensity distribution, hfr), may Be readily computed. 

Making use of the Fourier transform convolution theorem, we can 

write 


H(p) = GCp) (A-3) 

where capitalization of the functi£n 5 ...-indicates their zeroth order Hankel 
transforms, and the variabT‘e‘ ^ denotes spatial frequency. Explicitly, 
we have 


PCp) ~ ~~ 2 2 - • Iao ^ • exD 

+ Bog^J L II or^ P j 

+ Bog^ * exp [-Ti^Og^p^] I I 


and 


G(p) = Jq {2irap) 


(A-4) 


(A-5) 


where the symbol Jq denotes the zeroth order Bessel function. 

But, the radial intensity distribution of the image can be written 


hCr) = 2ir / p H(p) Jq (2irrp) dp 

Jq 


(A-6) 


Substituting equations (A-3), (A-4) and (A-5) into equation (A-6), we 
obtain 


h(r) = 




Tt + BOr 
I t 



(A-7) 


v.'here 



p Jq (ap) Jq * ^^‘P 


where 

a = 2ira 
3 = 27ir 


(A-8) 


(A-9) 


• 16 


Equation (A-8) is of a standard form tabulated by Gradshteyn and Ryzhik (1965). 
Explicitly, we have 


■^k = [“3 /^k J * /V3 * ^0 

°k 


where the symbol Iq denotes the modified Bessel function of the first kind 
of zeroth order. Formulae for the computation of are given by 
Abramowitz and Stegun (1972). 
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Figure 1 . 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


FIGURE CAPTIONS 

0 

Stellar Slit Scans (A6800A): Comparison of Calculation with 

Observation. Composite observational data are illustrated by 
the solid curve. The best-fitting analytical curve is illustrated 
by the broken curve. The slit width was 0.645 arc. sec. The 
sample interval was 0.032 arc, sec. 

0 

Stellar Slit Scans (X7300A): Comparison of Calculation with 

Observation. Composite observational data are illustrated 
by the solid curve. The best- fitting analytical curve is 
illustrated by the broker! curve. The slit width was 0.645 arc. 
sec. The sample interval was 0.032 arc. sec. 

Model Neptune Slit Scans. Theoretical Scans were computed for 
each of the six models listed in Table II. The Neptune disk was 
taken to be circularly symmetric with an apparent angular radius 
of 1.168 arc. sec. Atmospheric smearing was described by the' 
mean point spread function discussed in the text. The slit 
width was 0.645 arc. sec. All slit scans were normalized to 
equal total energy, obtained by integrating over the disk 
of the planet. Broadening and flattening of the slit scans 
increases progressively with the degree of limb-brightening 
on the disk. 

0 

Neptune Slit Scans (X6B00A): Comparison of Theory with Observa- 

tion. Composite observational data are illustrated by the 
individual dots. The sample interval was 0.032 arc. sec.. 
Theoretical scans, computed both for a uniform disk and for a 
ring of infinitesimal width, are illustrated by the solid 
curves. The Neptune disk was taken to be circularly symmetric 
with an apparent angular radius of 1.168 arc. sec. Atmospheric 
smearing was described by the relevant point spread function 
discussed in the text. The slit width was 0_^645 arc. sec. 

All slit scans, both observational and theoretical, were nor- 
malized to equal total energy, obtained by integrating over the 
disk of the planet. 

o 

Neptune Slit Scans (X7300A); Comparison of Theory with Observa- 
tion. Composite observational data are illustrated by the 
individual dots. The sample interval was 0.032 arc. sec. 
Theoretical scans, computed both for a uniform disk and for a 
ring of infinitesimal width, are illustrated by the solid 
curves. The Neptune disk was taken to be circularly symmetric 
with an apparent angular radius of 1.168 arc. sec. Atmospheric 
smearing was described by the relevant point spread function 
discussed in the text. The slit width was 0.645 arc. sec. All 
slit scans, both observational and theoretical, were normalized 
to equal total energy, obtained by integrating over the disk 
of the plant. 
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between two oppositely wound coils of the QRS inductor. To renormalize the 
oscillator's mass m one might attach to it a conducting plate that resides 
in the Inductor's magnetic field. The velocity of the plate through the 
magnetic field would indtice an electric dipole moment on the plate which, 
in turn, would couple by its velocity to the magnetic field giving on inter- 
action energy proportional to p n and thence a mass renormalization. 

Unfortunately, these various ideas have not yet produced a viable design 

for clean coupling of a mechanical oscillator's energy to a QRS. On the 

other hand, designs without clean coupling can still yield measurements of 

more accurate than the amplitude- and- phase quantum limit (n^ + •^)^ fiuJ* 

An example is a QRS that couples only to x , but that overages x over a number 

of cycles before sending it into the first classical si age (amplifier) (9, ll) . 

The measurement scheme of Fig, 1 will do this if the period of the circuit's 

(QRS) oscillations is much longer than the period of the mechanical oscillations. 

Then the circuit's capacitance (Eq. (24)) and resonant frequency will be 

sensitive to the time average of x and thence to H^, with only small admixtures 

of sensitivity to the time- varying part of x^ and thence to the oscillator's 

phase ([r. This is equivalent to the statement in Eq. (g) that Q^^ = 

f(fi + ct(H )ilf) with a « 1, which in turn permits accyraoiea much better than 
1 

AE = (n^ + |')2 ilui. Reference ( ll ) sketches a detailed analysis of this type 
of scheme, but for an electromagnetic oscillator with a mechanical QRS and 
with Rj. = Kx^q rather than Kx^q^ as in Fig. 1 and Eq. (24). That 
analysis reveals a limiting sensitivity 

1 X 

AE a (n^ + i)2 (n/u))2 hw, (25a) 

where E = (n^ + ^) Tnu is the oscillator's energy, (u is its frequency, fi is 
the frequency of the QRS, and n « w. The corresponding limit on the 
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detection of a classical force cos (u)t + cp), which drives changes In the 
oscillator's energy, is 


F 

o 


> 



(25b) 


if (i/n < n^ + •!■. If {i/n > n^ + 'J, then the limit (25a) on AE gets replaced 
by fw), the ultimate precision with which one can ever measure energy changes; 
and correspondingly die force limit (25b) gets replaced by (21). 

fsP A 

In measurements of ths time average of x' and thence it is not essential 

that the interaction Hamiltonian involve x^. Instead can be proportional 

^ 2 
to X, and then the internal workings of the QRS can produce the average of it 

at the entrance to the first classical stage. 


BACK-ACTION- EVADING MEASUREMENTS OF 

llie QNDF observable = x cos lot - (p/mw) sin lot (real part of complex 
amplitude; Eq. (5)), like the position x, has a continuous spectrum of 
eigenvalues; and in principle it can be measured arbitrarily quickly and 
accurately ( 13 . 15 ) . Suppose that an initial "state-preparation" measurement at 
t = 0 has put the oscillator into an eigenstate of Xj^(o) with eigenvalue 

5^. A classical force F(t) [total Hamiltonian H = H^ - x F(t)] will change 
X^ as seen in the Heisenberg Picture 

t 

X (t) = X.(0) - f [F(t')/nw)] sinuJt'dt' . (26) 

In the Heisenberg Picture the oscillator's state remains fixed in time at 
ll^), but this is an eigenstate of X^(t) with eigenvalue 

t 

|(t) ~ % - r [F(t')/moj] sinujt' dt' . 

° 0 


(27a) 



5 . 0 
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Subsequent perfect measurements of must yield this eigenvalue and will 
reveal the full details of its evolution. It evolves in exactly the same 
manner as would evolve for a classical oscillator ( 13 , 15) . 

One pays the price, in these meusurements, of not knowing anything about 
the imaginary part of the complex amplitude Xg (Eq„ (Sc)). However, if 
one has a second oscillator toupled to the same force F(t), one can mcanure 

/V 

the imaginary j^art Yg of its complex amplitude, giving up all information 
about the real part Yj^. One's measurements must give the eigenvalue 

t 

Ti(t) =115 + / [F(t*)/ms] cosout'dt' , (27b) 

which evolves in exactly the same manner as the Xg or Y2 of a classical 
oscillator. From the output of either oscillator, or better from the two 
outputs, one can deduce all details of the evolution of F(t), no matter how 
weak F(t) may be ( 15,15 ) . Thus, X^^ and Yg are QNDF observables. 

A perfect measurement of Xj^ (or Yg) requires (i) that the interaction 
Hamiltonian depend on X^^, and (ii) that commute with X^^ (Eq, (7) 
and associated discussion). The simplest example is 

Hj. = KXj^q = Kxq cos u)t - (x/muj) p q sinmt , (28) 

A coupling of this type can be achieved, for a mechanical oscillf.tor, 
by using a capacitive position transducer with sinusoidally modulated coupling 
constant = Kxq cos tut), followed by an inductive momentum transducer 
with modulated coupling constant (h^ = -(k/ihu) pq simot). The two trans- 
ducers together produce a voltage output 


9 = 3Hj/^q = Kx cos out - (k/cuu) p sin tut = KX^ , 


(29) 



5 . 0 
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which drives an electromagnetic circuit, the QRS, in which the charge q 
flows (15). While capacitive position transducers and inductive velocity 
transducers are easy to construct, inductive momentum transducers arc not. 

The momentum and velocity of the oscillator are related by 

X t= 5(H^ + H^)/cip = p/ro - (l^mu)) q sinwt, (30a) 

wliich means that the classicol Lagrangian L = px - (U^ + Ej.) for oscillator- 
plus- transducers is 

L = ^ mx^ - ^ mio^x^ - Kxq cos (Ut + (K/(u)xq ainuJt + -|r m(K sin u)t/nua)^q^. (30b) 

The first two terms represent the oscillator; the third is the capacitive 
position transducer; the fourth is an inductive velocity transduodr (wire, 
physically attached to oscillator, moves through external magnetic field); 
and the last Is a negative capacitor In the QRS circuit. Thus, an Inductive 
momentum transducer is equivalent to an inductive velocity transducer (easy to 
construct) plus a negative capacitor (hard) (15) , Although negative capacitors 

are not standard electronic components, they can be constructed in principle, 
and in principle they can be noise free (15) . 

For an electromagnetic oscillator with mechanical QRS, one can achieye 
the desired Ej. = K^j^q using a capacitive transducer for the oscillator's 

/s . . 

"position" X (= charge in oscillator s capacitor) and an inductive transducer 
for its "momentum" p (^ flux in oscillator's inductor). The momentum trans- 
ducer turns out to involve a standard mechanical current transducer (current s 
x) pl'is a negative spring in the QRS ( l5) . In principle negative springs can 
be noise free (15) . 

The sinusoidal modulations required in the transducers must be regulated 
by an external, classical clock, which has the same frequency lu as. one's 




S. 0 
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oscillator* One cannot use the oscillator Itself as one's clock because In 
e»;tractlng the required oscillatory Information from the osclllotor one wil 

A 

produce an unacceptably large back action on However, before the experi- 
ment begins one can check the frequency of one's clock against that of the 
oscillator. In principle they can be made to agree perfectly, one! In principle 
the clock can be made fully classical so Its outputs arc real numbers, cos uut 
and nlnujt, rather than operators ( lO, 15 ) . In pracMce, frequency drifts and 
quantum features of the clock need not cause serious experimental problems (15, 20) . 

A perfect meosurement of Xj^, which lasts a finite time t, requires in- 
finitely strong coupling in the transducers (K ->• oo) In order to give a signal 
that overwhelms zero-point noise in the QRS, If one has only finite coupling, 
then the zero-point noise accompanying the signal gives rise to a limit 
(Eqs* (16) and (17) with x ^ X^ and Nt ->■ t) (Refs, (13-15)) 

> (V2nuj)2 (Pujt)"^ . (31) 

Here P Is the dimensionless coupling constant (Eqs. (l6)). Thus, whereas 
stroboscopic measurements with limited coupling can beat the amplltude-and- 
phase limit by a factor of only (Pu)t)~^ (Eq« (I8b) with N = ojr/jt), continu- 
ous back-action-evading measurements of X^^ can beat it by (Pu)t)~^. Stroboscopic 
measurements are worse because of their smaller duty cycle. 

In the realistic case of weak coupling, P < 1, one must average over 
many cycles (ujt » l/ p) in order to substantially beat the amplitude and phase 
limit. In this case one can make use of a "trick" analogous to measuring the 
energy by coupling to x and averaging: One can perferm a "single- transducer, 

back-action-evading measurement" ( l3* l4, 15) by coupling to 


X cosmt cos fit s (Xj^ + Xj^ cos 2ujt + Xg sin 2ojt) cos fit 


(32) 
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(l.c*^ o 2Kx^ cos wt cos nt) and then sending the signal through a filter 
(the QRS) with bandpass at frequency n » u and bandwidth ^ « l/2% « u/2n. 
The filter will "average the Xg signal away" until its omplltudo has fallen 
by l/S(ST^ relative to that of the X^^ signal. Since the initial rms Xg signal 
strength is l/*/2 that of the X^, this corresponds to «= f (Xj^ + Xg/2^/2 u>T^.) 
in Eq. (6a), which together with the uncertainty relation (Sc) and the 
argument of Eqs* (6) tells us that (20) 

^ (h/SoKu)^ (2V2 iST.y.)'’^ . (35) 

This is the error in X, due to back action from measurement of X„, The ftddi- 

i, d 

tlonal error due to zero-point tolse accompanying the X^^ signal into the 
amplifier is (Eqs. (16) and (l7) with x -»• X^ and Nt t) (Refs, (l4). ( 20) ) 

fiXj^ > (VSnw)^' (^mT)■'2 . (:5k) 

Here l/Zx is the bandwidth of the experiment (t is the larger of the QRS 
averaging time snd the averaging time in subsequent electronics). Th»j 
ultimate quantum limit on the sensitivity is (33) if 0 > 2-j2x^x , and 
(34) if p < 2/2 t^7t. Note that Eq. (3*+) is the game limit (to 

within factors of order unity) as in the case of exact coupling to X^ = 

X cosiBt- (p/muj) sin tut. Thus, when {3 < 2*72 t^t and uir » 1, one can 
abandon the momentum transducer without any serious loss of accuracy! 

This type of single- transducer, time- averaged, back-action-evading 
measurement of X^^ appears today to be the most viable technique for beating 
the amplitude-and-phase limit (l4) in gravitational-wave detection. In place 
of (l4) one will face the limiting measurable force 

^ (^t) (muii) ^ • Max 



(35) 
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THERMAL NOISE IN THE OSCILUTOR AND AMPLIFIER 

The quantum Limits derived above arc not achievable in the laboratory 
today because thermal noise exceeds quantum mechanical noise. 

Ignore for the moment thermal (Nyquist) noise in the oscillator. Then 
if the resistors in the QRS arc cooled sufficiently, the dominant non-quantum 
noise will be that in the amplifier (first classical stage). The amplifier, 
which we assume to be linear, can be characterized by its power gain G and 
its noise temperature T^. The QRS feeds the amplifier a signal at frequency 
f = n/2it, to which the amplifier adds a noise power per unit bandwidth 

dP vr. 

n no 

df “ exp(i^kT^) - 1 ‘ ' ' 

Here k is Boltzmann’s constant. If the incoming signal has power P , then 

8 

the amplifiad signal and noise have power (25) 

° 3xV(r;| - ki„) 

Here Af is the bandwidth, and the hn/2 is a zero-point energy that accompanies 
the signal throughout its treck through the amplifier sni other electronics, 
but does not get omplified ( 25) . The quantum limits of previous sections of 
this article are attributable to this zero-point energy. In the presence 
of a real linear amplifier, with non-negligible noise temperature T^, the 
signal power P^ must fight not (fiD/2)Af, but rather 

( ^ X i Af 

[ exp(rjiykT^) - 1 G 2 ^ “ 

(cf. Eq. (57)), Consequently, it is reasonable to expect that the 
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omplifior nolae will modify our quantum limits (Eqs, (6c), (l3)-(lS), 
(17)-(19); (21), (22), (2S), (31), (33)-(35)) by replacing ft with (ig) 


« 5- 2ft ft 

^ n ^pTftn/kVy- 1 G 


2kT 


(38) 


n 


n 


if 


kT » ftn and G » 1. 
n 


These modified quantum limits are sometimes called "amplifier limits". 
Quantum mechanical analyses of linear amplifiers ( 2h ,25) reveal a minimum 
possible noise temperature (noise temperature cf "ideal linear amplifier") 


(T ) ^ 

n'min 



= ftn/k £n2 if G » 1, 


(39) 


which is achievable in principle by a maser amplifier, and which corresponds 
to 


Note that for large gains an ideal linear amplifier will give worse energy 
performance by a factor 2 (ft ->• 2ft) than the quantum limits of previous sec- 
tions. 

The best linear amplifiers that have been built are parametric amplifiers 
and maser amplifiers, which operate at microwave frequencies and have (kT^^^o) 
as small as ~ 10 ft. With such amplifiers one can only hope to get within a 
factor 10 or </l0 of our quantum limits (ft -> 10ft). And to even achieve this 
one must design a QRS which upconverts the oscillator's signal frequency 
(kilohertz in the gravitational-wave case) to the microwave (gigahertz) region. 
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Any physical oscillator (e„g., the fundamental mode of a gravitational- 
wave bar antenna) is weakly coupled to a thermal bath of dynamical systems 
(eig., sound waves in the bar). This coupling produces a frictional damping 
of large-amplitude motions, and it also produces a thermal-buffeting random 
walk of the oscillator’s amplitude ("Nyquist noise"). The root-raean-square 
random-walk change of the oscillator's amplitude during time t is 

^^o^Nyq "" ^^o/™^^Nyq “ ^^l^Nyq ~ ^^2^Nyq 

, pi „ , i 
~ (kT/mm )^ (tar/n)^ • 

Here T is the temperature of the thermal bath (the bar's temperature), and Q 
is the oscillator's quality factor (nurabcf of radians of oscillation required 
for frictional damping of large-amplitude oscillations by a factor e in energy). 
Che corresponding root-mean-square energy change is 

(AE)j.^q ~ (E^kT)2 (wt/q)^. (4 lb) 

These Nyquist noises must not exceed the amplifier limits (quantum limits with 
h -»■ 8kT^^^/fl) if one is to achieve the amplifier limits in real experiments. 

Some nuffibers w< 11 be given below. 

IT30SPECTS FOR STROBOSCOPIC MEASUREMENTS 

One possible scheme for stroboscopic measurements of a mechanical oscil- 
lator [gravitational-wave antenna with mass m 10 kg and frequency to ~ 

3 X 10^ sec is shown schematically in Fig. 2. The mass of the oscillator 
is physically attached to the central, movable plate of a capacitor [(capacitance 
between outer plates) = C], which plays the role of transducer. The capacitor 

1 ][Q 

resides in the QRS — a high-frequency LC circuit [frequency n = (LC,)~^ IS' 10 sec 
which has small losses [amplitude damping time t = 2(RCfl^) ^ « O.l/u)] and 
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which is driven at its resonant frequency fi by an external generator. In 
practice this circuit would be a microwave cavity; see Ref. (26). At the 
measurement times U)t s 0^ rt, 2it, ... the generator is turned on for a time 
f/2 and then turned off, — *,d in an additional time x/2 the excitations in 
the circuit die out. During the brief "on time" x, the amplifier sees a 
voltage signal V = (V /d) nxx cos nt, where V /d is the amplitude of the 
oscillating electric field between the cii.pacitor plates. The experimenter 
averages the amplitude of this signal (with alternating sign) over N 
measurements to determine the position x of the oscillator. 

It is straightforward to analyze the noise performance of this system 
using standard circuit theory. Alternatively, one can invoke the general 
formulas (15)-(19) for stroboscopic measurement schemes. Assuming that the 
resistor's physical temperature is less than "xplifier's effective 
temperature '’^eff - 10 K, the amplifier noise dominates and in Eq. (I8b) we 
must replace h -v 2kT^^^/n, Assuming that the amplifier is properly 
impedance-matched to the circuit, the measurement will achieve the limiting 
precision (I8b); 



Comparison of the voltage signal with Eq. (l6a) reveals that Klg^^^l = 

(V yj~2d) (Ht) f scrutiny of Fig. 2 reveals that the QRS output impedance, as 
seen by the amplifier, is gg^ = 2x/C; consequently the dimensionless coupling 
constant of Eq. (l6b) is p = (V ^d)^ Gfix/(^^mu)^). Combining this with the 

required pulse time x = (PNm ) ^ (Eq. (l8a)), we find 

PN = 
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To avoid voltage breakdowri in the capacitor, its electric-fieid amplitude 
should not exceed (V^d) ~ 10° volts/ cm. Assuming other reasonable parameters 
C 1 pf, n a 10^^ sec N m 1000, ~ 10 K, m 10 kg, uj rr 3 X 10^ sec 

we find 

PN cr 20, to ~ 1 X 10 cm. (44) 

Thus this system can achieve a sensitivity that is a factor (20)^rw2,l below 

the "amplitude-and-phase" amplifier llmitj} but this is still an order of magnitude 

worse than the amplltude-and-phase quantum limit (h/2moj)^ Ef I x 10 cm. 

Nyquist noise in the antenna (Eq. (4la) with lOr = xN) will be less 

-17 

than the measurement precision to ^ 1 X 10 cm if the antenna is 

Q 

cooled to 4 K and has a quality factor Q 2 ; 4 x 10" . This is comparable to the 
best mechanical Q that has been achieved ( 27) for a sapphire crystal at 4 K, 


PROSPECTS FOR SINGLE-TRANSDUCER BACK.- ACTION- EVADING MEASUREMENTS 

The configuration of Fig. 2 can also be used in a single- transducer, back- 

action evading measurement of In this case the circuit's amplitude damping 

2 -X 

time 2(RCn ) becomes the averaging time of the QRS filter (previously it was 
the stroboscopic pulse length); and we require t^, » l/w (previously it was 
« l/uj). Instead of being pulsed, the generator's modulating voltage has the 
steady-state form sinDt sin u)t, which produces an electric field 

(V ^ d) cos nt cos (jt in the capacitors (V^ = U^n/2u). That electi'ic field, 
interacting with the motions x = cos cot + siniot of the mechanical 
oscillator, produces a signal voltage 

Vg = (vyd)(nT^,./2) |X^ sin fit + (2uTj(.)'^Xj^ sin Dt sin2ut 

( 45 ) 

- ( 2 wt^.)”^X 2 sin fit cos 2u)t 
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at the output of the QRS. Amplification of this signal produces information 
about and Xg with relative accuracies £K^ = (2*>/2 * 

Assuming that the resistor noise ts negligible compared to amplifier 
noise (which it will be for the parameters used below and for a resistor 
physical temperature ^ ~ 10°K/^ we can compute the noise performance of 

this system ^rom Eqs, (33), (3lf), and (16) with li + SkT^^^t The best 

quality factor that hr.a been achieved ( 26 ) for a superconducting microwave 

n 

resonator (our QRS circuit) with a narrow capacitive gap is = OTj^. s; 10 , 
corresponding to ~ 10 sec. Consequently back-action forces (Eq. (33)) 

limit the sensitivity to 


AX, 


^ mw / 


(2/2 


= T 2 


X 10' 


■18 


cm. 


(46) 


a factor 7 below the amplitude-and-phase amplifier limit and appro?fimately 
twice the amplitude-and-phase quantum limit. (Here we use = 10°K, 

10-1 4-1 

£1 = 10 sec , ra = 10 kg, and lu = 3 x 10 sec as before.) In order that 

Nyquist noise in the mechanical oscillator (Eq^ (4la)) not exceed this 

sensitivity, the averaging time must not exceed r ~ 0.01 sec. (Here we ut-e 

o 0 

the same oscillator temperature and Q as before, T = 4 K and Q = 4 x 10 .) 

To achieve the limit (4 g) we also require a coupling constant 8 a? 2s/2/ai)T 2 ; 
0.005 (Eq, (34)). To compute Q, first derive Klgg^l = (V^d) (nT^/2«vf2) 
from 'Eqs, (4s) and (l6a) with x -v X^; then evaluate the impedance seen 
by the amplifier in Fig. 2 at the X^ signal frequency fi = (lc)"^^ - 2t^/C; 

then evaluate Eq, (l6b): 


8 = (V^d)^ cnTy(16mto^) . 


(47) 


The required 8 of 0.005 can be achieved with a reasonable electric field in 


Braginsky ~ 32 


the capacitive gap: V^d = 1 x 10^ volts/ cm. 

This example and that of the last section confirm that it is easier to 
achieve a given level of sensitivity by continuous, single-sensor back-action 
evasion than by stroboscopic techniques. However, along the route toward 
realization of such experiments tliere remain a series of difficult experi- 
mental problems ~ not least of which is the frequency stability of the clock 
that regulates the voltage generator. 

ON THE UMITING FREQUENCY STABILITY OF A GENERATOR 

Although current technology can achieve the frequency stability required 
by the above examples, it is of longer-term interest to know ultimate quantum 
mechanical limits on the stabilities of clocks. 

At present the world's most stable clocks are the "superconducting cavity 
stabilised oscillator" (SCSO) ( 28) and the hydrogen maser (29). Both involve se 
excited electromagnetic oscillations inside a cavity. In the SCSO the clock 
frequency (I is regulated by the cavity’s normal mode; and a change M of a 
typical dimension £ of the cavity will produce a frequency change 

An/n Ai/f. (4a) 

In the maser, if the electromagnetic quality factor of the cavity (teflon 
bubble) exceeds ^ 0 • (mean time hydrogen atoms spend in cavity) = Q^, then 
expression (40) will be true. Otheirwlae, AO/fi rr (Af/i) ‘ (Q^/Q^), and the limit 
derived below is correspondingly modified. 

Reference ( 50) derivea a qiiantum limit on the frequency stability of any 
electromagnetic oscillator satisfying (48). The source of the limit is 
quantum fluctuations in the deformation of the cavity walls by electromagnetic 
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fltressee. Since the b tresses in the electromagnetic fiold are equal to its 
energy density H^/ £ (with the Hamiltonian of the eiectroraagnetic oscil- 
lator), the force on the walls is and this deforms the walls by 5i « 

H /£k where k is the "mechanical spring constant" of the walls. The electro- 
magnetic field is in a thermalized coherent state with n^ quanta, which 
possesses quantum fluctuations s consequently, A£ J: n^ ^iin/k£, 

xfhich leads to frequency fluctuations (Eq. (48)) 






n. 


■hn/k£‘ 


(49) 


This "electromagnetic back-action" limit must be contrasted with the limiting 
precision for measurements of n during an averaging time t; An ^ where 

At ^ nj^ is the quantum uncertainty in the phase of the oscillator's coherent 
state; 

An/n n^“2 (nr)' (so) 


(" Townes -Schawlow limit"). These two limits lead to an optimal number of 
quanta n^, and an ultimate quantum limit 


n 





(r.i) 


For a cavity with wall thickness comparable to cavity dimensions £, or for 

a "cavity" made by coating the outside of a dielectric crystal with super 

conducting material ( 51) , the spring constant k is related to the Young's 

2 

modulus of the cavity walls by k — • v/ 1 , where V is the cavity voliune. 

Then 

An/n > (tt/E^VT)2. (51 «) 

In practice ^ 10^^ dynes/ cm, V — 1 cm^, so Afi/n Z 10~^^ • (x/l sec)"2. 

This limit is achievable in principle, but current technology is far from it. 
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Sunanary . Some future gravitational-wave antennas will be loO 
kilogram cylinders, whose end-to-end vibrations must be measured 
so accurately (l0“^^ centimeters) that they behave quantum 
mechanically. Moreover, the vibration amplitude must be measured 
over and over again without perturbing it ("quantum nonderaolition 
measurement"). This contrasts with quantum chemistry, quantum 
optics, or atomic, nuclear, and elementary particle physics where 
one usually makes measurements on an ensemble of identical 
objects, and one cares not whether any single object is perturbed or 
destroyed by the measurement. This article describes the new 
electronic techniques required for quantum nondemolition 
measurements, and the theory underlying them. Quantum non- 
demolition measurements may find application elsewhere in science 
and technology. 
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Fig. 1 


Fig. 2 


FIGURE CAPTIONS 


2 2 

« A scheme for coupling a mechanical oscillator's (position) ^ x 
to an electromagnetic QRS. 

. A scheme for stroboscopic or continuous back-action-evading 

measurements of a mechanical oscillator. This scheme was devised 
independently in 1978 by V. B. Braginsky and by R. W. P. Drever, but 
has not previously been published. 





